Coupling of actin hydrolysis and polymerization: Reduced descrip- 
tion with two nucleotide states 
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Abstract. - The polymerization of actin filaments is coupled to the hydrolysis of adenosine 
triphosphate (ATP), which involves both the cleavage of ATP and the release of inorganic phos- 
phate. We describe hydrolysis by a reduced two-state model with a cooperative cleavage mecha- 
nism, where the cleavage rate depends on the state of the neighboring actin protomer in a filament. 
We obtain theoretical predictions of experimentally accessible steady state quantities such as the 
size of the ATP-actin cap, the size distribution of ATP-actin islands, and the cleavage flux for 
cooperative cleavage mechanisms. 



Introduction. — Actin filaments are an important 
structural element of the cytoskeleton, and their ATP- 
driven polymerization dynamics plays an important role 
in cell motility [1] . The ATP hydrolysis in actin filaments 
is the basis for treadmilling, i.e., the simultaneous poly- 
merization and depolymerization at the two ends of a fila- 
ment, and is necessary for actin- mediated force generation 
and motility [2]. 

Actin monomers (G-actin) assemble into polar actin fil- 
aments (F-actin) with a fast polymerization dynamics at 
the barbed end and a slow polymerization dynamics at 
the pointed end. Actin monomers can bind ATP, which 
is then hydrolyzed in a two-step process into adenosine 
diphosphate (ADP) and inorganic phosphate (Pi). First, 
ATP is cleaved into the complex ADP-Pj, from which Pj 
is released in a second step. After incorporation into a 
filament the actin protomers can therefore be in three dif- 
ferent states: a T-state (ATP-actin), a 0-state (ADP-Pj- 
actin), and a D-state (ADP-actin). 

Past studies have focused on two sorts of cleavage pro- 
cesses in actin filaments. In random cleavage, T-protomers 
are cleaved independent of the state of their neighbors 
[SHS]. In vectorial cleavage, there is a sharp interface be- 
tween the ATP-cap containing only T-protomers and the 
remaining filament consisting of G- and D-protomers, and 
cleavage can only occur at the T0-interface [IHl]. 

In this Letter, we investigate cooperative hydrolysis 



mechanisms, where the cleavage rate of each monomer 
depends on the state of its neighbors and which contain 
random and vectorial mechanisms as special cases. Co- 
operative hydrolysis was previously discussed in Ref. [10] 
for actin and Ref. [Ill for microtubules. One important 
piece of evidence for a cooperative hydrolysis is the small 
hydrolysis rate of G-actin as compared to F-actin, which 
acts as a very effective ATPase with a fast hydrolysis rate. 
This pronounced change of the ATP hydrolysis rate af- 
ter inclusion of G-actin monomers into filaments suggests 
that the cleavage rate is affected by binding to other pro- 
tomers in the filament. Furthermore, structural differences 
between T- and D-state protomers indicate that cleavage 
might be affected by the same structural elements that are 
also involved in the binding of protomers [12[ll3j . 

A complete model of ATP hydrolysis involves all three 
nucleotide states of actin protomers. We have studied such 
three-state models for cooperative ATP hydrolysis in Ref. 
[131 . In this Letter, we will consider a reduced two-state 
model by combining the 0-state and the D-state of pro- 
tomers into a single D*-state. This means that we fo- 
cus on the cleavage of T-protomers and ignore the sub- 
sequent process of Pj-release. The reduction to two pro- 
tomer states can be justified for fast growth at high T- 
monomer concentrations. The reduced model has the ad- 
vantage that we are able to obtain analytic predictions for 
important steady state observables, such as the size of the 
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Fig. 1: Actin filament consisting of T-protomers (red) and D*- 
protomers (white) in the reduced two-state model, (a) Random 
cleavage with rate Uc; (b) vectorial cleavage with rate (jJc,d*', 
(c) cooperative cleavage with rate u>c,d* and uJc,t depending 
on the local neighborhood. 



ATP-actin cap, the size distribution of ATP-actin islands, 
and the cleavage flux. In particular, we can calculate ana- 
lytically how these quantities depend on the cooperativity 
of the cleavage mechanism. We find an intriguing scaling 
behavior in the limit of strongly cooperative cleavage. 

Tviro-state model for cooperative hydrolysis. 

In the reduced two-state model, we combine the 0-state 
and the D-state of protomers into a single D* -state and 
model the actin filament as a one-dimensional sequence of 
actin monomers in the T- or D*-state (ignoring the helical 
structure of the filament), see Fig. [T] 

In general, cleavage within the actin filament happens 
according to a cooperative mechanism, i.e., the cleavage 
rate of each monomer depends on the state of its neigh- 
bors. Each monomer is polar, therefore we do not expect 
cleavage and release rates to be mirror symmetric, i.e., 
they need not to be invariant under exchange of the two 
neighbors. The simplest way to introduce cooperativity 
without mirror symmetry is to assume that the cleavage 
rate depends on the state of the neighbor in the direction 
of one of the ends. Because the ATP-binding cleft is lo- 
cated in the direction of the pointed end, it is plausible to 
assume that cleavage depends on the state of the neighbor 
monomer on the "pointed side" . Consequently, we intro- 
duce two cleavage rates for an ATP monomer: Uct if the 
monomer has a T-neighbor on the pointed side and uJcD* if 
the monomer has a D*-neighbor on the pointed side. The 
two cleavage rates Wct, ^cD*, can also be written as 

UJcD* = Wc , WcT = t^cPc (1) 

which defines the cleavage parameter pc- We assume that 
the presence of a cleaved neighbor monomer increases the 
cleavage rate, which implies Pc < 1- In the special case 
of Pc = 1, the cleavage rate does not depend on the state 
of neighboring monomers, which corresponds to random 
cleavage. In the limiting case of pc = 0, cleavage only hap- 
pens at the TD*-interface within the filament correspond- 
ing to a vectorial cleavage mechanism. A small cleavage 
parameter pc ^ 1 corresponds to a strongly cooperative 
cleavage mechanism. 

In the following, we will focus on the barbed end of the 
filament. T-monomers attach with a rate Won (the number 
of T- monomers attaching per unit time) at the barbed end. 



Table 1: Literature values for model parameters. 
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This attachment rate is proportional to the concentration 
Ct of T-monomers in solution and the rate constant Kon, 
i.e., Won = KouCt- We assume that the attachment of 
D*-monomers is not possible, i.e., Won.D* = 0, which is 
justified in view of smaller rate constants |5| and small 
concentrations of Q- and D-monomers in solution. 

Both T- and D*-monomers can detach from the barbed 
end with rates Woff ,t and Woff .d* • Starting from the full 
three-state model, an effective detachment rate Woft.D* 
of D*-monomers can only be defined consistently for the 
whole T-monomer concentration range Ct if 6- and D- 
monomers have similar detachment rates, Woff^e — Woff,D- 
Measured values uJos^q — 0.2s^^ and uJos^b — 5.4s"^ [5] 
show that this condition is violated. Nevertheless, a re- 
duced two-state model can still be introduced for high T- 
monomer concentrations Ct, where the probability Pi^d* 
that the first protomer at the barbed end is a D*-protomer 
is negligible or the probability Pi^t that the first protomer 
at the barbed end is a T-protomer is close to one, 

Pi^D' ~ and Pi,T - 1 - Pi^D' ~ 1. (2) 

In this limit the corresponding detachment flux of D*- 
monomers -Pi,d* Woff,D* ~ is always negligible, and the 
detachment rates of both 0- and D-monomers become ir- 
relevant for the polymerization process. In the following 
we will focus on this limit Pi.d' ~ 0. We will show below, 
see eq. (fTTj) . that it is realized for strongly cooperative 
cleavage pc ^ 1 or for fast growth at high T-monomer 
concentrations Ct- 

We do not take into account a possible cooperativ- 
ity in the attachment and detachment process, i.e., Won, 
t^off.T, and Woff,D* do not depend on the state of the last 
monomer, which is at the tip before attachment or which 
is left behind at the tip after detachment. We also ne- 
glect fracture of filaments, which has been discussed for 
hemoglobin fibers in Ref. |16j . Literature values for cleav- 
age, attachment, and detachment rates of our model are 
listed in table [T] 

We will derive analytic results for the length of the 
ATP-actin cap, the length distribution of ATP-actin is- 
lands, and the cleavage flux. We compare these results to 
stochastic simulations of the full three-state model, which 
were performed using the Gillespie algorithm as described 
in Ref. |14j . For sufficiently high T-monomer concentra- 
tions Ct, we expect simulation results for the three-state 
model to agree with our analytic results for the reduced 
two-state model. In the stochastic simulations, we will use 
a three-state model with a random Pj-release mechanism 
with release rate ojr — O.OOSs"^ [3f8|. Furthermore, we use 
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^oB,e — 0.2s ^ [5 and WqS.d = 5.4s 
of O- and D-protomers, respectively. 



5 for the off-rates 



Growth rate and critical concentrations. At- 
tachment and detachment parameters for T-monomers 
lead to the T-protomer growth rate, 



Jt 



Pituj, 



off.T 



t^off'.T. 



(3) 



Because there is no D*-attachment, Won,D* = 0, and we 
consider the limit Pi_d* ~ 0, the T-protomer growth rate 
Jt equals the total growth rate Jg of the filament, 



Ja 



Pl,Tl^oS,T — -Pl,D''^off,D' ~ Jt- (4) 



The critical concentration for filament growth at the 
barbed end is given by Cr.g = i^off,T/'«on with Cr.g — 
0.12/^M for the values given in table [H 

For vectorial cleavage with pc — 0, there is a single 
ATP-island at the filament tip and a single TD*-interface, 
where cleavage takes place with a rate lucD* = <jJc- The 
length of this ATP-tip becomes infinite in the steady 
state if the T-protomers growth rate exceeds this cleav- 
age rate, Jt > which defines a threshold concentra- 
tion Ct,c = ('^c + '^off.T)/'*on with Ct,c — for 
the values given in table [TJ The dimensionless growth pa- 
rameter Jt/'-^c = (Ct — Cr.g) I {Ct,c — CT,g) characterizes 
the competition of growth and hydrolysis currents. Be- 
cause of condition ([2]), we will focus on fast growth with 
Jt/i^c ^ 1 in the following, see eq. ([TT|) below, which is 
realized for T-monomer concentrations much larger than 
the corresponding threshold concentration, Ct 3> Ct,c- 

Length distribution of ATP-tip. The growing 
barbed end consists of a sequence of T- and D*-monomers. 
The state of the filament can be described as a sequence of 
connected islands of T-monomers which are separated by 
D*-monomers. Similar to an analysis of the case of random 
hydrolysis in microtubules in Ref. [VT\ we will focus on the 
length distribution of these "ATP-islands" . First we will 
consider the length of the first ATP-island at the barbed 
end, which we call "ATP-tip" in the following. 

The probability pk of finding an ATP-tip of length k 
k = 0,1,2,..., with A; = corresponding to the case of a 
D*-monomer right at the barbed end, satisfies the master 
equation 

dtPk = JT(pfc-l -Pfc) ~ l^c [1 + /Oc(fc - - (^fc^o) 

+^cPk+l +^cPc ^ Ps (5) 
s>k+2 

with boundary condition p_i = 0. The first term on the 
rhs of eq. ^ describes loss and gain by attachment of T- 
protomers. The second term is nonzero only for tip lengths 
fc > and describes loss by cleavage: the last T-monomer 
at the TD*-interface is cleaved with a rate locD' , whereas 
the remaining fc — 1 T-monomers of the tip are cleaved 
with a rate uJct- The last two terms are the corresponding 
gain terms from cleavage: an ATP-tip length of k can be 
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Fig. 2: Length distribution pk of ATP-tip for an actin con- 
centration Ct ~ 1/iM (corresponding to Jt/loc ~ 34) and 
cleavage parameters pc = 10~^ (black, o) and 10""^ (red, □). 
Other parameter values as in table [1] Comparison between (i) 
analytic results from eq. ((8]) using pk = Pk — Pk+i (solid lines) 
and (ii) results from stochastic simulations using the Gillespie 
algorithm (circles). 



obtained by cleavage at the TD*-interface of an ATP-tip 
of length fc -|- 1 with rate lOcD* , or it can be obtained by 
"fragmentation" of a tip of length s > k + 2 into two pieces 
by cleavage of its interior T-monomers with rate uJcT- For 
the special case of random cleavage with pc — 1 we recover 
the results of Ref. pTT]; for vectorial cleavage with pc — 0, 
eq. ([5]) reduces to a random walk in /c-space (fc > 0) with 
stepping probability Jt from fc — 1 to fc and Wc from fc + 1 
to k. With a continuum approximation in the variable k, 
eq. (O has been obtained also in Ref. [Tl] . 

The steady state of the ATP-tip length distribution pk 
is obtained by setting the rhs of eq. ([5]) equal to zero. It 
is convenient to consider the cumulative quantity Pk = 
'^i>kPi '^ith ^0 = 1 a-iid Pi = Pi,T, for which we find the 
recursion relation 



k+i 



-Pk = 



Jt 



a;c(l - Pc) 



{Pk-Pk-i) + 



Pc 



1 



-kPk (6) 



in the steady state with fc > 1 . For vectorial cleavage with 
Pc — 0, we find an exponentially decaying stationary solu- 
tion Pk ^ {JtI^ci^ for Jt < t^c, i.e., below the threshold 
concentration Ct < Ct.c- For Jt > uJc, the ATP-tip is 
steadily growing and no stationary solution can be found. 

For the general case of arbitrary cooperativity, we apply 
a continuum approximation in the variable fc, 



= dlPkik) + adkPkik) - bkPk 

.1- Pc - Jt/^^c 



with 



a = 2- 



1 



Jt/oJc 



and b — 2pc 



1 



Jt/uJc 



(7) 



The solution of this equation with boundary condition 

Po = 1 is 



Pfc = cxp(— afc/2) 



Ai(b~2/3aV4 + 6i/3fc) 
Ai(6-2/3a2/4) 



(8) 



where Ai(a;) is the Airy function [18]. The ATP-tip length 
distribution is obtained as Pk = Pk ~ Pfc+i from the solu- 
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tion ([5]), see Fig. [51 The continuum approximation used 
to derive eq. ([7]) is justified because neither \a\ nor b can 
become large compared to unity. 

In the following we will focus again on fast growth with 
Jt/^c ^ 1, which leads to a « — 2 and b w 2pcUJc/Jt ^ 
1. Using the asymptotics of the Airy function, Ai(a;) ~ 
g-2x^/V3 for a; > 1 [ig, we find Pk ~ e-'''='/2|»l and, 
consequently, the resulting tip length distribution pk « 
—dkPk{k) is exponentially decaying for large fc, 

Pk « (fc6/|a|)e-'"'/'l'^l ~ k{Lu,p,/JT)e~^'^^'>^/^'-. (9) 

An interesting observable, which may be experimentally 
accessible in experiments on single filaments, is the mean 
tip size (fc) = X]fc>i kpk- For fast growth we obtain 
from eq. ^ a characteristic square-root dependence on 
the cleavage parameter pc, 

(fc) « ^{\a\/by^' « v^(^TKPc)'/^ (10) 

which could be used in experiments to determine pc by 
measuring the tip length. Another experimentally ac- 
cessible observable is the probability Pi^d* that the first 
monomer is in the £)*-state, for which we find 

Pl,D' = 1 - Pi « PcUJc/2Jt, (11) 

i.e., a linear dependence on pc- Eq. ([TT]) also confirms 
that the limit Pi,d* ~ 0, see ([2]), is attained for strongly 
cooperative cleavage pc ^ 1 or for fast growth Jt ^ cJc- 

Size distribution of ATP-islands. Experiments 
probing the structure of single filaments can give informa- 
tion not only on the length of the ATP-tip but the whole 
distribution of ATP-islands sizes in the filament. The av- 
erage number Ik of ATP-islands of length k fulfills the 
master equation 

dth = JriPk-i - Pk) - i^ci^ + {k - l)pc) h 

+Wc(l + Pc)Ik+l + iuJcPc ^ Is (12) 

s>k+2 

for fc > 1. The first term on the rhs of eq. ([72]) gives the 
change in ATP-island numbers from attachment and de- 
tachment at the tip. The second term describes the loss 
from cleavage at the fc — 1 sites of the ATP-island with 
T-neighbors on the pointed side with rate Wct and the 
loss from cleavage with rate tUcD* at the TD*-interface at 
the pointed side of the island. The third and fourth terms 
on the rhs of eq. (|12p are gain terms from cleavage at an 
interior site of an island. The fourth term means that any 
ATP-island of length k can be obtained by "fragmenta- 
tion" of an island of length s > fc + 2 in two ways with 
a rate lOcT for each way. The third term means that an 
ATP-island of length k can also be obtained by cleavage 
at the boundaries of an ATP-island of length k + 1 with 
rate uJcT on the barbed side and with rate uJcd-' on the 
pointed side of the island. 



Two important quantities, which follow from the ATP- 
island distribution and can be observed in experiments 
on single actin filaments, are the average total number of 
ATP-islands / = J2k>i and the average total number of 
ATP-actin protomers {Nt) = X]fe>i kik- The average to- 
tal number {Nt) of ATP-actin protomers gives a measure 
of the average total length of the ATP-cap of the actin 
filament. The average total number of ATP-islands / also 
gives the number of TD* -interfaces within the ATP-cap, 
where cleavage takes place with rate uJcD* ■ Knowledge of 
these two quantities therefore not only characterizes the 
filament structure but allows us to also calculate the re- 
sulting cleavage flux. For the two quantities / and {Nt), 
we obtain the rate equations 

dtl = JtPo + i^cPcUNT) - 21) - wc(l - Pc)Ii (13) 
dt{NT) = Jt - ujcPc{{Nt) - I) - ujj (14) 

Equation (fT3|) describes the creation of additional ATP- 
islands by addition of T-monomers at an empty tip or 
by cleavage in the interior of an existing T-island with 
rate WcT, whereas the island number is reduced by cleav- 
age of ATP-islands of unit length with a rate lOcD-' ■ The 
first term on the rhs of eq. ([T4l) describes the addition 
of T-monomers with the rate Jt, the last two terms the 
loss of T-protomers by cleavage: the total number of 
TD* -interfaces with a cleavage rate uJcD* is given by /, 
whereas the number of sites where cleavage happens with 
the slower rate uJct is given by {Nt) — I. Therefore, the 
total cleavage flux is given by 

J, = iJcPci{NT) - I) + ujj (15) 

and eq. ((Ti)) is equivalent to dt{NT) — Jt — Jc'- In a steady 
state we must have Jt = Jc, i-e., T-protomer addition flux 
and cleavage fiux balance. 

The steady state of the ATP-island distribution Ik ful- 
fills dtik — 0. We determine the stationary island distribu- 
tion for fast growth in two steps: (i) We obtain a differen- 
tial equation for Ik in k by requiring that dt{Ik^Ik-i) = 0. 
This will give the stationary Ik apart from one integration 
constant, (ii) We will determine this integration constant 
by looking for stationary solutions of cq. ([T4|) for {Nt). 

i) To determine the steady-state ATP-island distribu- 
tion from eq. we consider dt{Ik — Ik-i) for A: > 2 and 
use continuous fc, which leads to 

dtilk ~ 4-1 ) ~ - Jrdlpkik - 1) + (wc - i^cPc) dllk{k) 
- ujcPc{k ~ l)dklk{k) - iujcPJk- (16) 

The steady state fulfills dt{Ik ~ Ik-i) — 0. In the limit of 
fast growth, the term JTd^pk is exponentially small, and 
we find 

Ik = cil (v/2^(fc - 1)) . (17) 

The integration constant c/ has yet to be determined. 
The scaling function 1(2:) = 2~^^'^D_3{x/\/2)e^ where 
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Fig. 3: (a) Double-logarithmic ploat of the size distribution Ik of ATP-islands for an actin concentration Ct = ifiM (corre- 
sponding to JtI^c ~ 34) and cleavage parameters pc = 10~^ (black, o) and 10"'^ (red, □). Other parameter values are as in 
table [1] Comparison between (i) analytic results from eq. (|2ip (solid lines) and (ii) results from stochastic simulations using 
the Gillespie algorithm (circles), (b) Double-logarithmic plot of (A'^r) (blue), / (red), and I\ (black) in the steady state as a 
function of pc for Ct = {Jt/oJc = 34). Comparison between (i) analytic results from eq. (|2ip (solid lines), (ii) results from 

a stochastic simulations using the Gillespie algorithm (circles), and (iii) results from a numerical integration of the full master 
equations ([5} and (fT2|) (squares). 



D^{x) is Whittaker's parabolic cylinder function jT8], ful- 
fills the differential equation 



= 21" {x) - xl'{x) - 3X(x). 



(18) 



It decays as I{x) 



for x 3> 1, which gives rise to 



a power-law tail in the island distribution with the same 
scaling behavior 1^ ^ as for random cleavage [T71 . 
From the island distribution p7)) . we derive 



dxT{x) = 



ci 



{Nt}^ — / dxxl{3 

^Pc Jo 




(19) 



where we used the relations dxX{x) = — X'(0) = 1/4 
and dxxX{x) — 21(0) = a/tt/^j which follow from the 
differential equation (fT8|). 

ii) Using the results in eq. (ITU) and requiring sta- 
tionarity dt{NT) = to leading order in pc, we determine 
the integration constant 



This leads to the final results 



2 



Jt 



2 LUcy^ 



(20) 



(21) 



for the average total length (Nt) of the ATP-cap, the 
total number of ATP-islands / and the number /i of 
short ATP-islands of unit length, see Fig. [3l The result 
I Jrp/ujc ^ 1 for the total number of ATP-islands is 
remarkable because we have 1=1 for strictly vectorial 
cleavage. Thus a small cleavage parameter ^ con- 
stitutes a singular perturbation, which gives rise to a pro- 
nounced change in the ATP-cap structure with a dramatic 
increase in the total number of ATP-islands. 



ATP-cap length and cleavage flux. The scaling 
behavior of Ik in eq. (I?!]) leads to a characteristic scal- 
ing of island sizes k with the square root of the cleavage 
parameter pc, which determines the pc-dependence of ex- 
perimental observables such as the average total length of 
the ATP-cap length (Nt) and the cleavage flux Jc- 

According to the ATP-cap length (Nt) depends 
linearly on the T-protomer growth rate Jt and, thus, 
on the concentration Ct of T-monomers and increases 
(X 1/ for small cleavage parameters corresponding to 
strongly cooperative cleavage mechanisms, see Fig. IH^a). 
This result is corroborated by a scaling argument |14j . 
It shows that experiments on the total cap length in the 
limit of fast growth will allow to determine the cleavage 
parameter pc if the cleavage rate is known. On the other 
hand, the cap length in (j2ip only depends on the prod- 
uct l/ujc^/Pc — \J^c,D''^c,T- Thus, a mechanism with 
high cooperativity and large cleavage rate can give rise to 
a similar cap length as a random cleavage process with a 
low cleavage rate. This demonstrates that measurements 
of certain filament properties such as the ATP-monomer 
content of filaments do not allow to uniquely distinguish 
between a vectorial model with high cleavage rate [BHS] 
and a random model with lower cleavage rate [3H5]. 

From the coupled rate eqs. ([T3| and ([T4|) for / and (Nt), 
we can also calculate the characteristic time scale r ~ 
1/ujcy/fh- for relaxation to the steady state with stationary 
values (|211) . This experimentally relevant time scale is 

— 1/2 

increasing with the same pc -dependence as the steady 
state ATP-cap length for small cleavage parameters. 

According to eqs. (fT4| and (|15|) . the cleavage flux Jc 
and the T-protomer addition flux Jt must balance in the 
steady state, i.e., Jc = Jt- Because of eq. ^ the cleav- 
age flux depends then linearly on the T-monomer concen- 
tration Ct and becomes independent of the cooperativity 
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Fig. 4: (a) The total average ATP-cap length (A'r) and (b) the total cleavage flux Jc as a function of actin concentration Ct 
for different cleavage parameters Pc- Comparison between (i) analytic results for Ct 3> Ct.c (the critical concentration Ct,c 
is marked by an arrow) in the regime of fast growth (solid lines) and and (ii) results from a stochastic simulations using the 
Gillespie algorithm (data points). 



parameter for fast growth with Pi^t ~ 1, see Fig. [H^b). 
For vectorial cleavage, on the other hand, we have / = 1, 
and the cleavage flux is directly given by the cleavage rate, 
Jc = uJc and, thus, independent of the T-monomer concen- 
tration, see Fig. IH^b). This shows again that a non-zero 
cleavage parameter represents a singular perturbation of 
vectorial cleavage, and the cleavage flux is a sensitive quan- 
tity to differentiate between strictly vectorial cleavage with 
Pc ~ and strongly cooperative cleavage with small but 
nonzero pc- 

Conclusion. — In conclusion, we have introduced an 
effective two-state model for cooperative ATP-hydrolysis 
in actin filaments, where cooperativity is characterized 
by the cleavage parameter pc- The model contains ran- 
dom {pc = 1) and vectorial {pc — 0) cleavage as special 
cases. For this two-state model, we could obtain analytic 
steady-state results for quantities such as the size of the 
ATP-actin cap, the size distribution of ATP-actin islands, 
the total number of ATP-actin islands, which describe the 
structure of the actin filament, and kinetic quantities such 
as the cleavage flux. Measurements of these steady state 
quantities will allow to determine the cleavage rate and 
the cooperativity of the ATP-cleavage mechanism. Re- 
cently depolymerization experiments on individual actin 
filaments have become possible [5l[T9], which allow to ob- 
tain information on filament structure based on the differ- 
ent depolymerization rates of T-, Q- and D- monomers. 
The depolymerization experiments in Ref. |19j as ana- 
lyzed in Ref. [TJ and earlier kinetic data discussed in Refs. 
[71ITD] suggest that the cleavage parameter pc is as small as 
10~^ — 10~^ and that cleavage is thus strongly cooperative. 
Our analytic results become exact in the limit of small pc 
and fast growth Jt ^ Wc and, thus, can be directly ap- 
plied in this relevant parameter regime. Measurements of 
force- velocity relations for polymerization under force can 
also provide sensitive probes of the filament structure as 
discussed in Ref. 120 . 
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